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ABSTRACT. We consider affine systems in R" constructed from a given integral in- 
vertible and expansive matrix R, and a finite set B of translates, a^x := R~ 1 x + b; 
the corresponding measure fi on R ra is a probability measure and fixed by the selfsim- 
ilarity fi = |S| _1 X^bgB ^ oa b • There are two a priori candidates for an associated 
orthogonal harmonic analysis : (i) the existence of some subset A in R" such that 
the exponentials {e z ^' x }xeA form an orthogonal basis for L 2 (fi); and (ii) the exis- 
tence of a certain dual pair of representations of the C* -algebra Ojy where N is the 
cardinality of the set B. (For each N, the C*-algebra On is known to be simple; 
it is also called the Cuntz-algebra.) We show that, in the "typical" fractal case, the 
naive version (i) must be rejected; typically the orthogonal exponentials in L 2 (fi) 
fail to span a dense subspace. Instead we show that the C* -algebraic version of an 
orthogonal harmonic analysis, viz., (ii), is a natural substitute. It turns out that 
this version is still based on exponentials e l , but in a more indirect way. (See 
details in Section 5 below.) Our main result concerns the intrinsic geometric features 
of affine systems, based on R and B, such that fi has the C*-algebra property (ii). 
Specifically, we show that fi has an orthogonal harmonic analysis (in the sense (ii)) if 
the system (R, B) satisfies some specific symmetry conditions (which are geometric 
in nature). Our conditions for (ii) are stated in terms of two pieces of data: (a) 
a unitary generalized Hadamard-matrix , and (b) a certain system of lattices which 
must exist and, at the same time, be compatible with the Hadamard-matrix. A par- 
tial converse to this result is also given. Several examples are calculated, and a new 
w— ) , maximality condition for exponentials is identified. 

o 

1. Introduction 

The present paper continues work by the coauthors in [JP3-6], and it also pro- 
rS 

5— i 

vides detailed proofs of results announced in [JP4]. In addition we have new results 

going beyond those of the announcement [JP4]. We consider a new class of self- 
similarity fractals X, each X with associated fractal selfsimilar measure ft, such 
that L 2 (/u) has an orthogonal harmonic analysis in the sense of C*-algebras (see 
(ii) below). This possibility is characterized with geometric axioms on the pair 
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(X,fi). It is known since [St3— 4] that /U is typically singular (in the fractal case), 
and that in general only an asymptotic Plancherel type formula can be expected in 
the sense of [Bes]. Our present approach is based instead on C*-algebra theory. In 
particular, we use the C*-algebras O^ of Cuntz [Cu], and we give the orthogonal 
decompositions in terms of a dual pair of representations of On where N denotes 
the number of translations in the affine system which determines \i. 

For an orthogonal harmonic analysis, the following three possibilities appear a 
priori as natural candidates: 
(i) the existence of a subset A in MJ 1 such that the exponentials e\(x) := e lX ' x 

(indexed by A 6 A) form an orthogonal basis in L 2 (/i); 
(ii) the existence of a dual system of representations of some C*-algebra On say, 
(N = the cardinality of B), such that one representation is acting afnnely in 
x-space, and the other (dually) in frequency-space (where the frequency variable 
is represented by A in the above exponentials e ); and finally 
(iii) one might base the harmonic analysis on an orthogonal basis of polynomials 
in n variables obtained from the monomials x a := x^x^ 2 ■ ■ ■ x^ n (where a = 
(cki, . . . , a n ) is a multi-index, ccj = 0, 1, 2, . . . , 1 < i < n), by the familiar Gram- 
Schmidt algorithm. 

But it is immediate that both of the possibilities (i) and (iii) lack symmetry 
in the variables x and A. Moreover, it turns out that (i) must be ruled out also 
for a more serious reason. We show in Sections 6-7 below that, for the "typical" 
fractal measures //, none of the orthogonal sets {e\} in L 2 {ji) will in fact span 
a dense subspace. Specifically, there is a canonical maximally orthogonal {e\} 
system such that a finite set of "translates" (details in Section 5) of it does give us 
a dense subspace. It is this extra operation (i.e., "spreading out" the orthogonal 
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exponentials) which leads to our dual pair of representations of the algebra On- 

It also turns out that case (ii) is a natural extension of our orthogonality con- 
dition, studied earlier in [JP2] for L 2 (Q), now with S7 some subset in W 1 with 
finite positive Lebesgue measure, and L 2 (Q) considered as a Hilbert space with the 
restricted Lebesgue measure. For the case, when £1 is further assumed open and 
connected, we showed, in [JP2] and [Pe], that (i) holds (i.e., there is a set A such 
that {e^JAeA forms an orthogonal basis in L 2 {Vi)) iff the corresponding symmetric 
operators {\/— lg^-}? = i, defined on C^°(Q), have commuting selfadjoint operator 
extensions acting in L 2 (Q). It is well known that extension theory for symmetric 
operators is given by von Neumann's deficiency spaces. But, even when individual 
selfadjoint extensions exist for commuting symmetric operators, such extensions 
are typically non-commuting. Hence, we expect that, also for our L 2 ([i) analysis, 
there will be distinct symmetry conditions and selfadjointness conditions. 

For our present case, the pair (R, B) is specified as above, the affine maps are 
given by a^x = R~ l x + b, and indexed by points b in the finite set B. We get the 
measure fj,, and the Hilbert space L 2 (fi), by a general limit construction which we 
show must start with some L 2 (Q) example as discussed. But, for L 2 (p), we show 
that the analogous symmetry condition is related to a certain lattice configuration in 
MJ 1 (see Lemma 4.1 below), whereas the analogous selfadjointness now corresponds 
to a spectral pairing between B and a second subset L in W 1 , of same cardinality, 
such that the N by N matrix {e lb ' e }, (for b G B, £ £ L), forms a so-called unitary 
generalized Hadamard matrix, see [SY]. Then this matrix, together with the lattice 
configuration leads to a dual pair of representations, as sketched above and worked 
out in detail below. The two representations will act naturally on L 2 (fi) and provide 
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a non- commutative harmonic analysis with a completely new interpretation of the 

classical time- frequency duality (see e.g., [HR]), of multivariable Fourier series. 

When our "symmetry" condition is satisfied, we get a dual pair of self-similar 
measures, [Ib and Hl, and this pair is used in the proof of our structure theorem. 
Many examples are given illustrating when the "symmetry" holds and when it 
doesn't. A connection is made to classical spectral duality, see e.g., [JP1-3]. 

2. Basic Assumptions 

We consider affine operations in M. n where n is fixed; the case n = 1 is also 
included, and the results are non-trivial and interesting also then. A system s 
in W 1 will consist of a quadruple (R,B,L,K) where R € GL n (M), B and L are 
finite subsets in IR n , and both of them are assumed to contain the origin O in W 1 ; 
finally K is a lattice in M n , i.e., a free additive group with n generators. It will be 
convenient occasionally to identify a fixed lattice with a matrix whose columns are 
then taken to be a set of generators for the lattice in question. It is known that 
generators will always form a linear basis for the vector space (see e.g., [CS]); and 
it follows that the matrix is then in GL n (R). 

With the assumptions (to be specified), it turns out that we may apply Hutchin- 
son's theorem [Hu] to the affine system {(Jb}beB given by 

a b x := i? _1 x + b, x£l" (2.1) 

There is a unique probability measure \i on R n satisfying 



H=\B\~ 1 J2^oa^ 1 , (2.2) 

b€B 
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which amounts to the condition 

[ fdn = \B\- 1 J2 [ f°«bdfi (2.3) 

J b J 

for all jU-integrable functions / on M n . For the matrix R, we assume that some 
positive integral power of it has all eigenvalues in{AGC:|A|>l}, and we refer 
to this as the expansive property for R. (It is actually equivalent to the same 
condition for R itself.) The use of [Hu] requires the so called open-set- condition 
which turns out to hold when our system s has a symmetry property which we 
proceed to describe. We then also have the following compact subset X, defined as 
the closure (in R n ) of the set of vectors x with representation 



oc 
X = 

i=0 



Y^R-'h, b t eB. (2.4) 



If \B\ < | deti?|, where \B\ denotes the cardinality of B, then the fractal dimension 
of X will be less than the vector space dimension n of the ambient W 1 . (See e.g., 
[Ke] for details on this point.) In general, the measure fj, is supported by X, and 
we may identify L 2 (fi) with L 2 (X,n) as a Hilbert space. We will refer to X as 
the "fractal" even in the cases when its dimension may in fact be integral, and the 
"fractal" representation will be understood to be (2.4). Occasionally, we will write 
X(B) to stress the digit-set B. 

3. Generalized Hadamard Matrices 

The two sets B and L from the system came up in our previous work (see [JP2- 
4] and [JP6]) on multivariable spectral theory. The condition we wish to impose on 
two sets B, L amounts to demanding that the corresponding exponential matrix 

(e iM ) (3.1) 
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is generalized Hadamard, see [SY]. The term b ■ £ refers to the usual dot-product in 
W 1 . It will be convenient to abbreviate the matrix entries as, (b,£) := e l27Tb ' e . Since 
OgB and £ L by assumption, one column, and one row, in the matrix ({b,£))BL 
consists of a string of ones. Let the matrix be denoted by U: We say that it is 
generalized Hadamard if the two sets B and L have the same cardinality, N say, 
and if 

U*U = NI N . (3.2) 

It follows from this that then also UU* = Nljy. (This is just saying, of course, that 
the complex N by N matrix, N~ 1 ' 2 U is unitary in the usual sense.) 

We noted in [JP6] that the harmonic analysis of type (ii) is based on this kind of 
Hadamard matrices. (The matrices also have an independent life in combinatorics.) 
It turns out that the matrices are known for N up to N = 4. We will show, in 
Section 7 below, that this then leads to a classification of the simplest affine fractals 
(as specified) such that the analysis (ii) exists. We say that two matrices U of the 
form (3.2) are equivalent, if N is the same for the two matrices, and if one arises 
from the other by multiplication on the left, or right, with a permutation matrix, or 
with a unitary diagonal matrix. We now list below (without details) the inequivalent 
cases of type (3.2) for N < 4. (For higher N, such a classification is not known.) 
After our present preprint was circulated, we learned that the N < 4 classification 
had also been found independently, see references [Cr] and [Wer]. The purpose of 
our examples in Section 7 is to show how the equivalence classes of (3.2) lead to 
distinct examples of fractal measures fi, and how the different [/-matrices lead to 
different dual pairs of representations. 

We will postpone to a later paper a rigorous classification of the different systems 
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(R, B), and of the corresponding type (ii) harmonic analysis of L 2 (/i). But we feel 
that the N < 4 examples are sufficiently interesting in their own right. They also 
serve to illustrate the technical points in our (present) two main theorems. 

Notice the 2ir factor in the exponential (3.1) above. It is put in for technical 
convenience only. 

Remark 3.1. If we pick the string of ones as first row and first column, then the 
possibilities for U when N = 2 are 



1 1 
1 -1 



(3.3) 



for N = 3, 



1 1 1 

1 C C I (3-4) 



1 c c, 

where £ is a primitive 3rd root of 1; and for TV = 4, 



/ 1 1 1 1 

11-1-1 



1 —1 u —u 
\1 —1 — u u 



(3.5) 



where \u\ = 1, up to equivalence for generalized Hadamard matrices, see e.g. [SY]. 

4. Selfadjoint Systems 

Corresponding to the affine mappings (2.1) for a given system s = (R,B,L,K) 
we have 

n (t) -=R*t + £, iGf" (4.1) 

and the inverses 

T7\t):=R*~\t-£) (4.2) 
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where the translations for (4.2) are given by the vectors 

b' ■= -R*' 1 (£), as £ varies over L. (4.3) 

Here R is an n by n matrix as specified above, B and L are finite subsets in W 1 
both containing O, and K is a rank n lattice. The invariance R(K) C K will be 
assumed, and we summarize this by the notation K 6 lat(i?). 

We introduce the dual system s° defined by s° = (R* , B' , L' , K°) where K° is 
the dual latttice, 

B':=-R*~\L) (4.4) 

and 

L' := -R{B) (4.5) 

The system s is said to be symmetric if 

R{B) C K, (4.6) 

and if K G lat(-R); and it is said to be selfadjoint if both s and s° are symmetric. 
(Also notice that, in general, we have s°° = s when s is an arbitrary system.) 

(The definitions are analogous to familiar ones for closed operators S with dense 
domain in Hilbert space, see e.g., [Fu]: the operator S is said to be symmetric if 
S <Z S* , where S* denotes the adjoint, and the inclusion refers to inclusion of graphs. 
It follows that S is selfadjoint, i.e., S = S* , iff both S and S* are symmetric.) 

We shall need the fact that B embeds into the of coset space R _1 (K) / K when 
additional orthogonality is assumed: 

Lemma 4.1. Consider a system s = (R,B,L,K) in M. n with the matrix R, the 
two finite subsets B and L in W 1 , and a lattice K as described above. Assume that 
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L C K° and that the two conditions (3.2), i.e., that Hadamard property, and (4.6) 
hold, where K° is the dual lattice in W 1 . Then it follows that different points in B 
represent distinct elements in the finite group R~ x (K) / K . 

Proof. Suppose b ^ b' in B. Then ^ e&L {£,b - b'} = using (3.2). But, for all 
k £ K, we also have X^ei(^j^) = l-^l> an d it follows that b — b' £ K; i.e., the 
R~ 1 (K)/K cosets are distinct. □ 

We shall assume in the following that our given system is of Hadamard type, 
i.e., that \B\ = \L\ and that the matrix (3.1) formed from (B,L) is generalized 
Hadamard, see (3.2) above. 

The following lemma is also simple but useful. 

Lemma 4.2. A given system s = (R, B, L, K) in W 1 is selfadjoint if and only if 
the following three conditions hold: 

(i) K G lat(-R), 
(ii) R(B) c K, and 
(iii) Lc K°. 

Proof. A calculation shows that K G lat(i?) holds iff K° E lat(i?*). For the system 
s° to be symmetric, we need R*(B') C K°, and that is equivalent to (iii) by virtue 
of formula (4.4). So both s and s° are symmetric precisely when (i)-(iii) hold. □ 

Remark 4.3. (Classical Systems) In [JP2], we considered the following spectral 
problem for measurable subsets Q C M n of finite positive Lebesgue measure, i.e., 
< m(tt) < oo where m = m n denotes the W 1 - Lebesgue- measure: Let Vt be given, 
when is there a subset A C 1" s.t. the exponentials 

e x {x) = (\,x)=e^ x - x , (4.7) 
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indexed by A € A, form an orthonormal basis in L 2 (fl) with inner product 



m(ft)- 1 / f(x)g(x)dx? (4.8) 



The problem (in its classical form) goes back to [Fu], and it is motivated by a 
corresponding one for commuting vector fields on manifolds with boundary, see also 
[Jol-2], [Pe], and [JP2]. 

We showed that the general problem may be "reduced" (by elimination of "triv- 
ial" systems) to a special case when the pair (Q, A) is such that the polar 

A° = {tGM": (t, A) = 1, VA 6 A} (4.9) 

is a lattice in R n , say K := A°, and the natural torus mapping W 1 —> W 1 / K is then 

i-i on n. 

In this case, there is a system s = (R, B, L, K) which is self-adjoint and of 
Hadamard type. Moreover the set A (called the spectrum) may be taken as 

k = L + R*K°. (4.10) 

Pairs (ft, A) with the basis-property are called spectral pairs; the "reduced" ones 
where A may be brought into the form (4.10) (with L ^ {0}, i.e., \L\ > 1) are called 
simple factors. We showed in [JP6] that more general ones may be built up from 
the simple factors. 

The following easy fact will be used below: Let K\ and K 2 be lattices, and let 
K i and K 2 be corresponding matrices. Then we have the lattice inclusion K\ C K 2 
if and only if the matrices factor: K\ = K 2 M with M E Mat n (Z) where Mat n (Z) 
denotes the ring of integral n by n matrices, i.e., M = (mij)™ i=1 with rriij G Z. 
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This observation allows us to take advantage of the Noetherian property of the ring 
Mat n (Z). A minimal choice for K subject to conditions is then always well defined. 
For a given lattice K, the dual lattice is denoted K° and given by 

K° :={s£R n : s-k£Z, V/c G K} 

If K is a matrix for K, then the inverse transpose, i.e., (K tr )~ l will be a matrix 
for K°. 

When R G GL n (M) is given, we denote by \&t(R) the set of all lattices K in W 1 
such that R(K) C K. For the matrices, that reads K _1 RK G Mat n (Z). This fact 
will be used in the paper; it implies for example that | det R\ is the index of K in 
R- X (K). It is known (see e.g., [CS] or [JP6]) that, if lat(-R) ^ 0, then deti? G Z. 



(Remark: If R is not in diagonal form, i.e., 



fr ... 0\ 

r ... 

\0 ... r/ 



rl n for r G Z, then 



there are lattices K not in lat(i?).) 

The standing assumption which is placed on R is referred to as the expansive 
property: We assume that, for some p £ N, all the eigenvalues {Aj} of R p satisfy 
\Xj\ > 1. Recall, R has real entries, but the eigenvalues may be complex. For 
emphasis, we will denote the transpose of R by R* , even though it is the same as 
R tr . (Note that the assumption on the eigenvalues of RP for some positive power p 
is equivalent to the same condition on R itself, i.e., to the condition for p = 1.) 

5. Iteration Systems 

In this paper, we shall study fractals (in the sense of (2.4) above) with a high 
degree of symmetry; and show that these fractals are precisely those which may be 
built from systems s = (R,B,L,K) which are selfadjoint, of Hadamard-type, and 
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where the lattice K is chosen as minimal relative to the three conditions (i)-(iii) 
in Lemma 4.1. In describing our limit systems (typically fractals), we show again 
that the Hadamard condition (3.2) is the central one. 

Motivated by (4.10), we form the set C(L) consisting of all (finite) sums 

£ + R*£ 1 +R* 2 £ 2 + --- + R* m £ m (5.1) 

when m varies over {0, 1,2,...} and £i G L. Using (4.1), also notice that C{L) is 
made from iterations 

n {nA---{n m m ■■■))■ (5.2) 

The set A in (4.10) is \J{Te(K°) : £ G L}. We shall also need the corresponding 
iterations, 

UKo...or c (F):! ie i}. (5.3) 

TO 

For a given string (£ , • • • , £ m )i the set in (5.3) will be denoted K°(£ 0l . . . , £ m ). 

Definition 5.1. We say that K° formed from a given system s = (R,B,L,K) is 
total if the functions {e s : s G K } C L 2 (fi) span a subspace which is dense in the 
Hilbert space L 2 (fi) defined from the Hutchinson measure fi, see (2.3). 

Both of our main results will have the total property for K° as an assumption. 
The way to test it in applications is to rely on our earlier paper [JP2] about spectral 
pairs, i.e., subsets £1, and A, in W 1 such that £1 has finite positive n-dimensional 
Lebesgue measure, and the exponentials {eA : A G A} form an orthogonal basis for 
L 2 (Q). We show that, for every such pair, the set 

K := A° = {i E R n : £ • A G Z, VA G A} 

is a lattice. Analogously to the situation in Lemma 4.1 above, we also show in 
[JP2] that the set fi in a spectral pair embeds in the torus W 1 /K. We identify a 
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special class of spectral pairs, called simple factors which produce two finite sets 
B, L C K n , and a matrix R with K G lat(R) such that the system s = (R, B, L, K) 
satisfies the conditions from section 4 above. Our present paper is motivated by 
getting "invariants" for simple factors from iteration of the affine maps (see (2.1) 
and (4.1) above). In [JP6] we further study the converse problem of reconstructing 
simple factors from "fractal" iteration limit-objects. In any case, the fractal limit 
X(B) from (2.4) will also be embedded in the torus W 1 / K. When equipped with 
Haar-measure L 2 (K n / K) has the exponentials {eA : A G K } as an orthogonal 
basis. In testing for our totality condition relative to L (//), we can then use that 
X(B) is the support of fj,, and then apply Stone- Weierstrass to {e\}\ € x° when 
viewed as a subset of C{X(B)). 

We shall say that C{L) is maximal if {eA : A G C{L)} is orthogonal in L 2 (/j,) and 
(considering t G W 1 ) if 

whenever (e t , e\) M = jl(\ — t) = for all A G C{L), then t G C{L). (5-4) 

We have used the transform jl given by 

fi(s) = I e s dfi = I e i2ns - x dfi(x) for s£R n . (5.5) 



We say that the system s is A-orthogonal, if the functions {eA : A G A} are or- 
thogonal in L 2 (ji), here A is given by (4.10), i.e., A = L + R*K° . (See also (5.6) 
below.) 

We are now ready for the 

Theorem 5.2. Let s = (R,B,L,K) be a selfadjoint system in M n , and assume 
(i) K° is total; 
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(ii) C(L) is maximal in L 2 ([i); and 

(iii) s is A- orthogonal, i.e., the points in A from (4-10) are orthogonal for £ ^ £' 
in A. 

Then it follows that s is of Hadamard type; i.e., \B\ = \L\ and the B / L-matrix U 
satisfies (3.2). 

Proof. Condition (i) states that the orthogonal complement of {e s : s € K°} in 
L 2 (ji) is zero. Notice that condition (iii) is equivalent to: 

fi{£ - £' + R*s) = 0, W^'inL, VseT. (5.6) 

If we set 

B(t):=|B|- 1 5^(6,t), VteR n ; (5.7) 

beB 

then (2.4) implies the factorization: 

jl(t) = B(t)fi(R*~ 1 t). (5.8) 

For distinct points £ and £' in L we claim that 

R*-\£-£')£K°. 

Assuming the contrary, there would be some s G K° such that £ — £' = R*s. From 
the orthogonality property (5.6) (see Definition 5.1), we then get 

A(0) = fi{£ -£'- R*s) = {e e , e e , +R « s )^ = 

contradicting /}(0) = 1. 

Since C(L) C K°, the point t := R* {£ — £') is not in C{L). From the maximal 
property (5.4), we conclude that there is some A E C(L) such that 
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This term works out to 

/i(ir _1 (*-o-A)(^o). 

From the orthogonality (Definition 5.1), we also have: 

= jX{£ - £' - R*X) = B{£ - 0£(-R* -1 (£ ~ O ~ A) 

where the last factor is non-zero. It then follows that B(£ — (!) = 0. 

Recall, for u,v € M n , the notation (u,v) := e l2 ' KU " v . Then the vectors {{■,£)} are 
indexed by points £ £ L, and we showed that they are orthogonal when viewed as 
elements in £ 2 (B). 

It follows that \L\ < \B\ where the symbol | • | denotes cardinality. We claim that 
they are equal. For suppose the contrary, viz., \L\ < \B\. Then pick coefficients 
kb £ C, not-all zero, indexed by b £ B, such that 

J2h(b,£)=0, V£eL. (5.9) 

b£B 

For every s £ K° and £ £ L, consider t := £ + R* s; and define 

f := Yl foXVb+R-HX)) ( 5 - 10 ) 

b£B 

where x denotes "indicator function" , the subscript is a 6-translate, and finally X is 
the i?-fractal. (Recall, details below, it is compact, and satisfies X = B + R~ l (X), 
with 

^((6 + J R- 1 (A > ))n(fe / + J R- 1 (X))) = for all b ^ b' in B.) (5.11) 

Note, (5.11) is a consequence of the totality of K° and the following observation: 
if b and c are in B and b + R~ l x = c + i? _1 y, then R(b — c) = y — x; it now follows 
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from Lemma 4.1 that x € y + K. But then 



(f,e t )v = / f(x)e t (x)dfj,(x) 

= £*>/ «<» + ,) 4*,) 

y^hf (£ + R*s,b + x)dfi(x) 

t~i Jr-Hx) 



b€B 



= [y>6<*,&>] / Mm (5.12) 

N v ' 

(5.9) 

(where we use (5.11) and Lemma 5.3 below), and 

(R*s,b) = (s,Rb) = 1. 

The last fact is from axiom (4.6) which makes Rb G K. It follows (from (5.9)) that 
/ is in the orthogonal complement of {ee-\-R* s } as £ varies over L, and s over K° . 
But from (i), we know that this is a total set of vectors in L 2 (fj,), so the function 
/ must vanish identically, [i-a.e. If the coefficients {kb} are not all zero, this would 
contradict (5.11), (2.3), and the basic properties of the Hutchinson measure \x. 

From the contradiction, we conclude that \L\ = \B\; which is to say, both con- 
ditions on the matrix ({b,£)) b », indexed by B x L, to be of generalized Hadamard 
type, are satisfied. We have \L\ = \B\ = N. If the matrix is denoted U, then 

UU* = U*U = NI N (5.13) 

where In denotes the identity matrix in iV variables, and U* is the transpose 
conjugate. To define it, it is convenient to use a common index labeling, e.g., 
{1,2,. ..,N}. □ 



HARMONIC ANALYSIS OF FRACTAL MEASURES 17 

Remark. Note that if we assume B is a subset of a set of representations for 
R~ l K J K , then (5.11) follows from an application of [Ke, Theorem 10] and a re- 
lated result in [Ma]. (See also [Ba-Gr] for related work.) The Kenyon-Madych 
result applies in the present context since the mapping from the set of all finite B- 
strings (bi, . . . ,b m ) with m varying in N, b t € B, into ^2 i R l bi G K is 1-1. This 
follows by induction and use of our orthogonality assumptions. To use [Ke]-[Ma], 
we then extend B so as to get a full set of residue classes R~ l (K) / K . 

Question. Does either of the following two conditions imply the other: (i) K° is 
total, (ii) B is a subset of a set of representatives for the quotient (R~ 1 K)/K1 

In the calculation (5.12) above, the following lemma was used. (It is needed 
because we do not know if, in general, [i is a Hausdorff-measure.) 

Lemma 5.3. Under the assumptions of Theorem 5.2, it follows that 



f(x) dfi(x) = f(x + b) dn(x) 

a b X JR- 1 X 



for all b in B and f in L 2 (/j,). 

Proof. The claim is equivalent to having 

fi(a b A) = M(i? _1 A) 

for all /^-measurable sets A C M n and all b in B; which in turn is equivalent to 

/i(a b A) = n(a c A) (5.14) 

for all ^-measurable A, and all b and c in B. The last equivalence used the assump- 
tion that € B. By regularity, it suffices to consider the case where A is a closed 
set. 
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Let A be a closed subset of X(B) = X, and choose B k C B k = B x • • • x B (k 
terms) such that 

(&!)■■■ i bk,b k+1 ) C B k+1 => (b 2 , ■ ■ ■ , b k+1 ) 6 5 fc , 



and such that 



A = f| E k , 



k=\ 

where 

E k = [J cr bl • • -a bk X. 

(6i,... ,b k )€B k 

The first condition means that (-E^t) is a descreasing sequence of compact sets. 
It follows (analogously to (5.11)) that the overlaps in the definition of E k are /i- 
null-sets. Hence, to prove the lemma, it suffices to show that 

M(cr&i ■■■a bk X) = n(a Cl ■ ■ ■ a Ck X) 

for all (pi,... ,b k ), (ci,... ,c k ) in B k . 

First note that, for any Borel set A, and any b in B, we have 

MA) = IBI" 1 E^c _1 A) > l^rVK'A), 

and hence ^(cr b A) > |i?| _1 ^(A). From this inequality, and (5.11), it follows further 
that 

H(X) = //( [J a 6 X) = J] n{a h X) > M (X), 

feG-B beB 

and therefore that ^{a^X) = |S| _1 /i(X). Assuming 

M ((7 62 ---C7 6fc X) = | J B|- fc +V(^), 
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it follows (analogously to the above), that 

A»(<7 6l <7 6a • • • °b h X) = \B\- k n(X). (5.15) 

Hence, by induction, (5.15) is true for all positive integers k, and all 6i, . . . , bk in 
B. This completes the proof of the lemma. □ 

Note also (heuristically) that (5.14) is a consequence of (5.11), and that 
^(o-fcA) = lim V] \B\~ 1 x<j b &{crb 1 ---Ob k x) 

bi,... ,b k 



lim S^ \B\ 1 XA((&b 2 ■■■<Tb k x) + Rb 1 - Rb) 

k—>oo *—* 



bi,... ,b k 

where x in X is arbitrary. However, the first equality requires that Xa b A is contin- 
uous. We refer to [Fa, p. 121] for further details on this point. 

We conclude with the following lemma which is both basic and general; in fact it 
holds in a context which is more general than where we need it. Such more general 
contexts occur, e.g., in [St4], [Mat], [MOW], and [Od], (among other places). But 
we will still restrict the setting presently to where it is needed below for our proof 
of Theorem 6.1. 

Lemma 5.4. Let (R, B) be an affine system in M. n (see details in Section 2) with 
R expansive and B C W 1 a finite subset. Let B be given by (5.7), and let \x be 
the 'probability measure from (2.2). We are assuming the property (5.11). Let 
J\f := {t G R n : B(t) = 0}. Then, for the roots of jX, we have 



{t 6 R n : flit) = 0} = \jR* k (M). 



fc=0 



Proof. We have (5.8) by virtue of [JP6, Lemma 3.4], and it follows that fi(t) = 

k 

when t E R* (M) for some k £ {0, 1, . . . }. From the assumed expansivity of R, we 
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also know that the corresponding infinite product formula is convergent, see [JP6, 

fc 

(3.13)]. In fact lim^ _> oo jl(R* t) = 1, for all t £ W 1 . This is from continuity of fi, 

k 

and the limit, R* t — ► 0. Now consider, 

fc-i 

3=0 

— k 

and suppose /}(£) = 0. Pick k (sufficiently large) s.t., fi(R* t) ^ 0. (This is 
possible by continuity, and the fact that A(0) = 1). We conclude, then that, for 
some j, < j < k, R* t £ A/"; and this is the assertion of the lemma. □ 

6. Orthogonal Exponentials 

We keep the standing assumptions on the quadruple s = (R,B,L,K) which 
determine a system in M n . In particular, the matrix R is assumed expansive (see 
section 2), the sets B and L in M. n are finite both containing 0. We will assume 
now that s is selfadjoint and of Hadamard type. We say that the system s is 
irreducible if there is no proper linear subspace V C M. n (i.e., of smaller dimension) 
which contains the set B, and which is invariant under R, i.e., Rv £ V for all 
v £ V. If such a proper subspace does exist, we say that s is reducible. In that 
case, it is immediate that the fractal X{B) from (2.4) is then contained in V . All 
the examples in Section 7 below can easily be checked to be irreducible. But the 
following example in M 2 is reducible, and serves to illustrate the last conclusion 

from our theorem in the present section: Let R = ( 9)'^ = i(n)'(n 
L=l(°),( 1 \\,K = Z 2 -,andV= ( ^ j , i.e., the x-axis in JR 2 . Then a 
direct calculation shows that X(B) = I I , where / = [0, 1] is the unit-interval on 
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the x-axis; and the Hutchinson measure [i is the product measure Ai % So where 
Ai is the restriction to / of the one-dimensional Lebesgue measure, and 5q is the 
point-O-Dirac measure in the second coordinate. For the set C(L) from (5.1), we 
have, 

£ ( l >={U«>) : " =0 ' 1,2 ""}' 

n represented by finite sums, n = YIjxj^Cj, £j £ {0, 1}, and N{n) = Ylj>o J^~ le j- 
We will show below that if K° is total then in the "fractal case", i.e., when 
N = \B\ = \L\ is smaller than |deti?|, and when s is irreducible, then C(L) is 
maximally orthogonal, and s is A-orthogonal. Recall [i is the Hutchinson measure, 
see (2.3), C{L) is the set given by (5.1); and finally the properties regarding the 
two sets K° and C(L) refer to the corresponding exponentials e s , when s is in the 
respective sets, and each 

e s (x) = (s,x) = e l27TS - x (6.1) 

is considered a vector (alias function on W 1 ) in the Hilbert space L 2 (/j,) = L 2 (X, /i). 
We refer to Definition 5.1 and Theorem 5.2 for further details. Recall that the 
total property (iii) in Theorem 5.2 amounts to the A-orthogonality, including the 
assertion 

/!{£-£' + R*s) = (6.2) 

for all I ^ I' in L and all s 6 K° . But if s is selfadjoint and of Hadamard type, 
then (6.2) follows immediately from (5.8), which is the functional equation of the 
transform fx, see also (5.5). 

The purpose of the present section is twofold. First we show that Theorem 
5.2 has a partial converse, and secondly that the technical conditions from our 
two theorems 5.2 and 6.1 amount to the dual pair condition (see Section 1) for 
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representations of the C*- algebra On- This is for systems s = (R,B,L,K) as 
specified where the two given finite sets B and L in R n are assumed to have the 
same cardinality N, i.e., \B\ = \L\ = N. Our recent paper [JP6] further details 
how the representation duality relates to our present assumptions. But we shall 
summarize the essentials here for the convenience of the reader. The Cuntz-algebra 
On (see [Cu]) is known to be given universally on N generators {s^ and subject 
only to the relations: 



s*Sj = Sijl and /_^ s i s *i = 1 (6-3) 



This means that, if a finite set of N operators Si say, acting on some Hilbert space 
TC say, are known to satisfy the relations (6.3), then there is a unique represention 
p of On, acting by bounded operators on TC, such that p(si) = Si for all i; or, 
equivalently, p(a)f = af for all a € On and all / € TC, where the operator a is 
given by the same expression in the SjS as a is in the Sj-generators. 

For a given system s = (R,B,L,K), there is then the possibility of making a 
representation duality based on the exponentials e l27Tt ' x in (6.1), and treating the 
two vector- variables x and t symmetrically: The pair (R, B) gives one affine system 
cjftX = R _1 x + b (b €E B) in the x-variable; and the dual system (R*,L) given by 
T(t = R*t + £ (I 6 L); a second one now acting in the t-variable. See (2.1) and (4.1) 
above. To be able to generate the asserted representation pair we need to specify 
{rt}e£L for enough values of t such that the corresponding functions e t (x) := e' l27Tt ' x 
span a dense subspace in L 2 {p). But it turns out that other conditions must 
be met as well: For the two affine systems {(Jb}beB and {t(}^l, the question is 
if we can associate operator systems {<S&} and {Z^} of 2N operators acting on 
L 2 (p,), each system satisfying (6.3), and the operators collectively defined from the 
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exponentials e t27Tt ' x as specified. When this is so, we have an (orthogonal) dual pair 
of representations of On acting on L 2 (fj>), and conversely. Then it turns out that, 
for each pair (b,£), the operator S^Tg is a multiplication operator, see (6.6) below. 

Theorem 6.1. Let s = (R,B,L,K) be a system in R n and assume that s is 
selfadjoint and of Hadamard type. Assume further that K° is total (with a minimal 
choice for K), and that \B\ < | det R\, and let n be the corresponding measure (see 
(2.3)) with support X. Then s is A-orthogonal and carries a dual pair of Cuntz 
representations (with On acting on L 2 (ji) for both representations) . If s is also 
irreducible, then C(L) is maximally orthogonal. 

Proof. Let N = \B\ = \L\ and note that from [JP6] (Theorem 4.1) we get a dual 
pair of representations {Sb}beB and {T^i^l of the Cuntz algebra On, see also [Cu] 
and [Ar], acting on L 2 (fj,) and given by the respective formulas: 

S * f = N -l/2 fo(Jb for /€L 2 (/U)) beB (64) 

T£e s = e Ti ( s } for s G K° and £ 6 L. (6-5) 

Moreover S^T^ is the multiplication operator Mm on L 2 (ji) given by 

M M f = N- 1 / 2 (e i oa b )f for / G L 2 ^). (6.6) 

It follows from [JP6, Theorem 4.1] that for £ and £' in L 

T e *T e ,=5u>I and ^TiTf = I (6.7) 

Here we use the Kronecker delta notation 



,1 ]££ = £' 

1 o Mi^e, 



24 PALLE E.T. JORGENSEN AND STEEN PEDERSEN 

and / denotes the identity operator in the Hilbert space L 2 (/j,). 

It follows then from (6.6) that the vectors e^_|_^* s , s £ K°, are mutually orthog- 
onal in L 2 (fi) for distinct values of £, i.e., for £ ^ £' in L. For more details on this 
point, we refer to sections 3-4 in [JP6]. For 

n 

A = ^iT%, (6.8) 

3=0 

we have 

Te T tl ---Tt n eo = ex. (6.9) 

Let A = Y^j=o-R* Iji an£ i K = YlT=o R* kj) where the £j's and fej's are in L. 
Then e\ and e K are orthongonal in L 2 (p) except in the cases where m < n, and 
£j = for j > m, and where m > n and fcj = for j > n. In the exceptional 
cases, it follows from (6.9) and T^e® = eo that e\ = e K . To prove the orthogonality 
assertion above, note that 

(eA,e K ) M = (Te T ei ■ ■ ■ Te n e ,T ko T kl ■■■T km e ) 

is = 0, unless £q = k , because T k Tj» = if k ^ £ . If k = £ , then T)* Tn = I, 
and we can repeat the argument on £\ and k\. It remains to consider the case where 
n > and m = 0; in this case, we will use the identity, T eo = eo, to write 

(e\, e K ) M = {T £o T ei ■ ■ ■ Ti n e ,T e )^ = 

and we conclude that (e\, e K ) M = unless £o = £\ = ■ ■ ■ = £ n = 0. 
It follows that the map, 



(£ ,... ,£ n )^Y, R * J ^ eC (^ ( 6 - 10 ) 
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is 1-1 on the set of finite sequences (£o, . . . , £ n -i,£ n ) with n a nonnegative integer, 
the ij's in L, and £ n ^ 0. 

The assumption that s be irreducible is now imposed, and we show that C(L) 
has the stated maximality property: We show that, if t G W 1 and {e\,e t )^ = 
for all A € £(L), then it follows that t G C(L). We shall do this by contradiction, 
assuming the t fi £(L). We shall use the functional equation (5.8) for p,, recalling 
that 

(ex,e t ) IJl =(i(t-\). (6.11) 

We shall also use that for every s G K n there is some £ <E L such that B{£ — s) ^ 0. 
This follows from the formula (5.7) for /?(■), and from the Hadamard property (3.2) 
which is now assumed. 

As a special case of (5.8), we get 



= fi(t - £ - R*h) = B(t - £ )KR* 't-R* X £ Q - £ t ) 



Picking £ G L s.t. B(t - £ ) ^ 0, we get 

= ^(R^t - R*^£ -h- R*£ 2 ) 



— i , — i 



= B{R* t-R* £ -h)p,(R* t-R* £ -R* £i-£ 2 )- 

Picking £\ G L s.t. 

B{R*~' t - R*~' £ - £i) ^0, 

we conclude next that 

fi{R*~ 2 t - R*~ 2 £ - R*^£i - £ 2 ) = 0, 
and we continue by induction, determining £ ,£i, . . . G L such that the points 



s p := R* P £ + ■■■ + R*^£ p -i + £ p 
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are in the dual fractal set X(L), see (6.2) above. When N < | deti?|, we may pick, 

inductively, the "digits" £i such that the differences 



,-p 



R* t-s p (6.12) 

are distinct as p varies, but 

fi(R*~ P t-s p ) = and B(i2*" ( "" 1) t - s p -i) + 0. 
Notice that the analytically extended transform 

fr(z) = I e l2nz ' x dfj,(x) (6.13) 



is entire analytic on C n , where for z = (z±, . . . , z n ) E C n , z-x = z\X\ + • • - + z n x n is 
the usual dot-product. Hence its zeros cannot accumulate. But the "dual attractor" 
X{L) (see (7.2)) is compact in W 1 so there a subsequence s Pi with limit s Pi — > s E 
X(L), and 



= lim/i(i?* l (t) - s Pi ) = A(-s) = /i(s) 



p- 



contradicting that the roots of /}(■) must be isolated (see (6.12)), even isolated in 
C n . The contradiction completes the proof, and we conclude that C{L) is maximal. 

If only a finite number of the "digits" £j are nonzero, then, using the contractive 

— 1 —p 

property of R* , we see that the sequences R* (t), and s p , both converge to zero 
as p — > oo, contradicting that /}(0) = 1, since A — ► p,(X) is continuous on W 1 . 

Claim 1. The set B is a subset of a set of representatives for R _1 (K) / K . 

Proof of Claim. From the self-adjointness of s we have RB C K (by Lemma 4.2). 
Therefore, -B C R~ l K. If 6 and b' are distinct and both in B, and if 6 E 6' + -ftT, 
then e t (b + x) = e t (b' + x) for all x E i? -1 ^ (all t E if ), contradicting the totality 
ofF inL 2 (/i). 
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Claim 2. The finite set L is a subset of a set of representatives for K° / R*K° . 

Proof of Claim. By Lemma 4.2, L C K°. If £ and £' are in L, and £ = £' + R*j for 

some 7 6 if°, then 

(b,£) = (b,£' + R*i) 

= (b,£')(Rb, 7 ) 
= (bj') 

where the last equality used Lemma 4.2 again. But this contradicts the Hadamard- 

property, unless 7 = 0. Considering, 

x p = R*~ P t - (R*~ P £ H + R^'ip-i + £ p ) , and 

y p = R* P x p = t-(£ + R*£i + ■■■ + R* P £ P ) , 

and letting P = {p : £ p ^ 0}; then we showed above that P is infinite, and that 
p € P — > y p is a 1-1 map. Hence {y p : p £ P} is infinite. □ 

Remark 6.2. For the reducible example (in K 2 ) mentioned in the beginning of the 
present section, we note that all the conditions of the first part of Theorem 6.1 
are satisfied. We also described the set C{L) of orthogonal exponentials for the 
example. But the maximality condition is not satisfied relative to L 2 (fi). Indeed, 
for the transform fi(s) from (5.5), we have, with s = (si, S2) £R 2 , 

| e ^!!£M if ( Sl ^0) 

A*(« =i -, Sl * - f V ^ „ 

[ 1 11 si = 0. 

It follows that the identity from (5.4) will be satisfied whenever £ = (£1, £2) £ H^ 2 
is such that £1 G Z_, i.e., negative and integral. (Specifically, /t(A — £) = for 
VA € C(L).) From the calculation of C(L), we note that such points £ = (£1,^2) 
will not be in the set C(L); and so the maximality condition is not satisfied for the 
example. 
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7. Examples 

7.1 Background Material. We now give examples to illustrate the conditions in 
Theorems 5.2 and 6.1. Since the generalized Hadamard matrices are known up to 
N = 4, the examples we give are "typical" for the possibilities when N < 4, and it is 
likely that there is a classification; but as it is unclear what is the "correct" notion 
of equivalence for systems s = (R, B, L, K) we will postpone the classification issue 
to a later paper. Note that the examples occur in pairs, one for s and a dual one 
for s° . Also note that each s will correspond to a spectral pair (£1, A) as well as 
a selfsimilar iteration limit, typically a "fractal" X with a self similar measure \i. 
When the given system s is selfadjoint, then there will in fact be a pair of "fractals" 
occurring as iteration limits, a selfsimilar jjl from the affine system: 

a b x = R~ l x + b leading to X = X(B) (7.1) 

defined from s, and also 

T f 1 (t) = lV~\t-i) leading to X (L), (7.2) 

and defining the corresponding dual selfsimilar measure //. Recall both jjl and // 
are probability measures on W 1 ; \i is determined by (2.2), and fj,' by: 



n 1 = \L\- 1 Y^^'°n, (7.3) 



see also (7.1)-(7.2) and Lemma 4.1 for more details on the dual pair of affine 
systems. 

Our examples below will be constructed from the matrices (3.3)-(3.5) which we 
listed in section 3. In fact, we shall supply a group of examples for each of the 
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generalized Hadamard matrices N = 2, N = 3, and N = 4, all the examples will be 
symmetric and of Hadamard type; but some will not be selfadjoint. In fact, when 
considering s = (R,B,L,K) we shall fix the first three R, B, and L, but allow 
variations in the lattice. When we insist on the Hadamard type, we shall see that, 
in some familiar fractal- examples, it will then not be possible to choose any lattice 
K such that the corresponding system s = s(—,K) is selfadjoint. We will then say 
that the system is not self-adjoint; it turns out that the obstruction is a certain 
integrality condition; and, when it is not possible to find a lattice consistent with 
both selfadjointness and Hadamard type, then it will typically be a simple, case by 
case computation, and we shall be very brief with detailed calculations. (It will be 
immediate that each of the examples in the list is irreducible; see section 6.) 

7.2 Group 1 Examples. We take N = 2; the matrix is (3.3), and the examples 
are illustrated with subsets of the line, i.e., n = 1, for W 1 . First, take R = 4, i.e., 
multiplication by the integer 4; the sets B and L will be B = {0, 1/2}, L = {0, 1}, 
and lattice K = Z. The (fi, A) spectral pair will be as follows: 

Q = [0, 1/4] U [1/2, 3/4], (i.e., the union of two intervals); 

A = {0, 1} + 4Z, (i.e., two residue systems modulo 4, see (4.10) above for the 
general case) 

X = iteration fractal, see Figure A, fractal dimension D = j^-| = 1/2, see the 
affine system (2.2), and also more details on fj, in Section 2 of [JP6]. 

It is easy to check that with this choice for R, B, L, and K, the corresponding 
system s is selfadjoint and of the Hadamard type. For this particular example, 
there are only two choices for K such that the corresponding system sk = s(—, K) 
is selfadjoint. They are K = % and K = 2Z. But the following modification, 
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corresponding to the classical middle-third-Cantor set, will only be symmetric; not 

selfadjoint: With 

R = 3, 

B = {0,2/3}, and 

L = {0,3/4}, 

we have the Hadamard type, c.f., (3.2); but there is no lattice K in R which makes 
the corresponding system sk selfadjoint (Graphic illustration, Figures A and 1). 

7.3 Group 2 Examples. We take N = 3; the matrix is (3.4), and the examples 

are illustrated with subsets of the plane R 2 . Take 
6 N 



R '0 6 I ' 



-'MMi- 



Z-M°W where /= § f 1 



if = 3Z 2 , i.e., multiples of the unit-lattice in 2 dimensions, equivalently 



/3m N 

points in R of the form I I where m, n £ 

3n 



The corresponding system will be selfadjoint of Hadamard type. If K is instead 
taken to be the lattice generated by the two vectors ( x ) and (3) (which turns out 
to yield K C L°), then there is a corresponding spectral pair (ft, A) where Q is 
a suitable union of scaled squares in the plane, and the corresponding spectrum 
satisfies A = K. But with this K, the iteration system s K will not have K° 
total in L 2 (/j,). In all, there are only three distinct choices, in this case, for lattices 
K in R 2 such that the corresponding system sk = s(— ,K) is selfadjoint: They 
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K\ C K2 C K3 for the lattices. The fractal dimension of X is 



In 3 



(Graphic illustrations, Figures 2-9.) 

7.4 Group 3 Examples. We take N = 4; the matrix is (3.5) corresponding to 
u = — 1, and the examples are illustrated with solid sets, i.e., pictures in 3-space 

R 3 . Take 

2 0' 
Z? - I 2 

2 

0\ /-|\ / \ / 



D ^ I 1,1 , I -I |.| 



K = l? = K° (i.e., selfduality). 

It is convenient to summarize the choices for R, B, L and K as follows: 

R = 2/3, 



S = -^h, and 



1 1 

III 1 

1 1 



where 



1 
1 
1 





-il 


' 1 


1 


-1 


L~ l = 


v-1 


-1 

1 


1 

1 
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is the unit-matrix. The choices for the lattice K are subjected to the conditions 
in Lemma 4.1. It turns out that the choice K = H? is the minimal one such that 
the system sk = s(—,K) is selfadjoint; and there is also a unique maximal choice 
for K with Sk selfadjoint, viz., K = L° where L° is given by (7.4) below. (Since 
L is symmetric, the matrix for the lattice L° is L~ x .) The lattice L° has matrix 
represented by the inverse 



(7.4) 



The corresponding system s = (R, B, L, K) is selfadjoint of Hadamard type. 
If the choice for K = 1? is replaced by K = L°, then (the modified) s is still 
selfadjoint: Notice that K = H? is the minimal choice for K (subject to (i)-(iii) 
in Lemma 4.1); and K = L° is the maximal one. This means that K° = 1? is 
maximal among the possible choices for K°; and this K° is total, see Definition 
5.1. 

7.5 Dual Pairs. The fractal dimension is D = |^-| = 2 which is integral, but less 
than the dimension (viz., 3) of ambient IR 3 . The fractal for the system s arises from 
scaling iteration of the set f2 = union of 4 cubes, see the figure (Figure B). For 
the dual system, Q° is instead the union of tetrahedra resulting in a 3-dimensional 
Sierpinski gasket, same fractal dimension D = 2, but with angles 60° rather than 
90°. The sketch is Figures 10-17, see also [Sch] for more details; it is the Eiffel 
tower construction, (maximal strength with least use of iron.) 

The corresponding planar Sierpinski-gasket corresponding to R = I I , B = 

\lii and L = 3 f I , does not have a lattice choice for K which makes the 
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associated system s = (R, B, L, K) in M 2 selfadjoint. The fractal dimension is 
D = M wl.58. 

In 2 

For the matrix (3.5) with primitive 4th roots of 1 (e.g., u = i), there is also a 
realization in R 3 : We may take R = 2I 3 , B = |l 3 , and L = 1 1 will give 

VI i \) 

a system of Hadamard type in 3-space, but again there is no lattice choice for K 
in M 3 such that the corresponding sk is selfadjoint. (Graphic illustrations, Figures 
10-17.) 



8. Concluding Remarks 

The operators {T#} from (6.5) and (6.7) may also be used in the definition of an 
endomorphism 9 on a certain C*-algebraic ©Ar-crossed product, it say. It is given 
by, 

6(A) = Y^ T eAT e * , for A e il , 

and clearly, 0(A*) = 0(A)*, and 0(AB) = O(A)0(B) for all A,B G it. Contin- 
uous versions, also called endomorphism-semigroups, have been studied recently 
by Arveson and Powers, see e.g., [Ar]. As spectral-invariants for these, Arveson 
has proposed (in [Ar]) a Cuntz-algebra construction which is based on Wiener- 
Hopf techniques, and which is inherently continuous, in fact with R + used as index 
for the generators in place of the usual finite (or infinite) discrete labeling set 
{1, . . . ,A^}. For our present B/L duality project with dual fractals, X(B) and 
X(L); we plan (in a sequel paper) to study an analogous C*- algebra construction 
which is generated by X(L) in place of IR+, but still modelled on Arveson's Wiener- 
Hopf approach. It appears that such an A(L)-fractal-based C*-algebra will serve 
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as a spectral-invariant for our B/L Hadamard-systems which are only symmetric, 
but generally not selfadjoint (relative to some choice of lattice K, see Section 7 
above) . 

The spectral-invariant question is an important one, and in our case we produce 
the dual representation pair (6.4) and (6.5) as a candidate. But representations {Sb} 
of On m the form (6.4), without a paired dual representation {T^}, cf. (6.6), are 
present for iteration systems which are much more general than the affine fractals 
studied here. As a case in point we mention Matsumoto's [Mat] recent analysis of 
(von Neumann type) cellular automata (details in [MOW] and [Od]); it is based on 
an ^-representation which is given by a formula similar to our (6.4) above. There 
is also an associated endomorphism with an entropy that can be computed; but we 
stress that for these (and many other) iteration systems, there is typically not a 
dualitly based on exponentials e and typically not a second {T#}- representation 
such that the two form a dual pair in any natural way. 

We have studied the class of spectral systems s = (R, B, L, K) in W 1 with special 
view to the selfadjoint ones which are also of Hadamard type, see Lemma 4.1. 
(When s is given in this class, the two sets B and L then have the same cardinality; it 
will be denoted N for convenience in the following comments.) It is important (but 
elementary) that this class of systems is closed under the tensor- product operation; 
i.e., if si and s 2 are systems in M ni and W 12 respectively, then the two properties 
(selfadjointness and Hadamard type) carry over to the system si <S> s 2 hi M™ 1 " 1 "" -2 = 
R ni x W 12 . 

If the Hilbert spaces for the respective systems are L 2 (/ii), i = 1,2; then the 
Hilbert space for si (g> s 2 is L 2 (/ii (g> /i 2 ), and the measure H\ (g> [i-i is the unique 
probability measure on IR" -1 x W 12 which scales the affine tensor operations of S\%s<i, 
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see (2.3) above. The set B for si <S> s 2 is B\ x B2, and the matrix-operation is, 
(61,62) 1— * (R\bi, i?2&2)- I n verifying the Hadamard property (3.2) for Si (8) $2; we 
use the important (known) fact that the class of generalized Hadamard matrices 
is closed under the tensor-product operation, i.e., that U\ ® U-i satisfies (3.2) with 
order N = NiN 2 if the individual factors Ui, i = 1,2, do with respective orders iVj, 
i = l,2. 

We say that a system s is irreducible if it does not factor "non-trivially" s — 
s\ <8> S2; and we note that the examples above from Section 7 are all irreducible in 
this sense. (In fact this irreducibility notion is different from that of Section 6, but 
the examples are irreducible in both senses.) 

The spectral geometry for regions in R n has a long history, see e.g., [Bo-Gu], 
[CV], [Ge], and [Gu-St]. But, so far, the Laplace operator has played a favored 
role despite the known incompleteness for the correspondence between the geome- 
try of the given domain and the spectrum of the corresponding Laplace operator. 
The approach in [De] is based instead on a multitude of second order differential 
operators, but the spectral correspondence is still incomplete there. Our present 
approach leads to a complete spectral picture and is based instead on a system of 
first order operators. For the fractal case however, the differential operators have 
no analogue. 

While our simultaneous eigenfunctions are based, at the outset, on a commu- 
tative operator system, our spectral invariant derives instead from a dual pair of 
representations of a certain non-abelian (in fact simple) C*-algebra. 

Self-similar limit constructions have received much recent attention, starting with 
[Hu], and then more recently, see e.g., [Ba-Gr], [Ed], [Ma], and [Ke]. These results 
seem to stress the geometry and the combinatorics of the infinite limits, and not 
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the spectral theory. Our present emphasis is a direct spectral/geometry- correspon- 
dence; and we also do not in [JP6] impose the strict expansivity assumption (which 
has, so far, been standard almost everywhere in the literature). Furthermore, we 
wish to stress that the sets £1 C M. n which occur in our present spectral pairs are 
more general than the self-reproducing tiles (SRT) which were characterized in [Ke, 
Theorem 10]. However, Kenyon's SRT's can be shown to satisfy our conditions, al- 
though our class is properly larger; not only because of the expansivity assumption, 
but also because of the combinatorics, see [Jo-Pe5] for details. Further work on 
these interconnections is also in progress. 
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